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We report inelastic neutron scattering experiments performed to investigate the low energy mag- 
netic excitations on single crystals of the heavy-fermion superconductor PrOs4Sbi2. The observed 
excitation clearly softens at a wave vector Q = (1,0,0), which is the same as the modulation vector 
of the field-induced antiferro-quadrupolar ordering, and its intensity at Q — (1,0,0) is smaller than 
that around the zone center. This result directly evidences that this excitonic behavior is derived 
mainly from nonmagnetic quadrupolar interactions. Furthermore, the narrowing of the linewidths 
of the excitations below the superconducting transition temperature indicates the close connection 
between the superconductivity and the excitons. 

PACS numbers: 74.70.Tx, 78.70.Nx, 74.70.Dd 
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The relation between magnetism and superconductiv- 
ity in heavy-fermion (HF) systems is one of the most 
interesting research areas in condensed matter physics. 
Magnetic fluctuations due to strong Coulomb repulsion 
and hybridization between f electrons and conduction 
electrons are believed to be the origin of both the HF be- 
havior and the superconductivity. In general, however, 
f electrons in systems with high crystal symmetry also 
have nonmagnetic orbital degrees of freedom and they 
can interact with conduction electrons through the or- 
bital channels. Therefore, the nonmagnetic orbital fluc- 
tuations might also lead to the HF superconductivity. 

The filled skutterudite compounds RT4X12 (R: rare 
earth or actinide, T: transition metal, X: pnictogen) crys- 
tallize in the body-centered cubic space group Jm3, in 
which R ion is surrounded by a cage of 12 X ions. These 
compounds show a wide variety of thermal, magnetic and 
transport properties due to strong c—f mixing effects and 
a 'rattling' motion of R ion resulting from the unique 
crystal structure. Among them, PrOs4Sbi2 is the first 
Pr-based HF superconductor with the superconducting 
transition temperature T c = 1.85 K Q, Q. The large 
Sommcrfcld coefficient 7 ~ 350 mJ/mol K 2 of the specific 
heat and the large jump at T c (AC/T C ~ 500 mJ/mol K 2 ) 
suggest that heavy quasiparticles participate in the su- 
perconducting transition. Specific heat, magnetization 
and thermal conductivity show double superconducting 
phase transitions No coherence peak in 1/Ti has 

been found at T c in the NQR experiment Further- 
more, ^SR measurement reveals that the superconduct- 
ing state breaks the time reversal symmetry |6j . The exis- 
tence of a field-induced antiferro-quadrupolar (AFQ) or- 
dered phase above ~5T near the superconducting phase 
naturally suggests that the nonmagnetic quadrupolar de- 



grees of freedom play an important role in this unusual 
superconductivity |2lla|- In PrOs4Sbi2, 4f electrons are 
well-localized under the cubic crystal field (CF) with 
Th symmetry. This is evidenced by dHvA experiment 
which shows the closeness of the Fermi surface between 
PrOs4Sbi2 and LaOs4Sbi2 and by inelastic neutron 
scattering (INS) experiments which show the clear CF 
excitations Q, 0, ^| . Compared to other HF supercon- 
ductors, the well-localized nature of 4f electrons is char- 
acteristic to this material. The CF level scheme is well es- 
tablished by elastic Q and INS experiments 0, |n| and 

(2) 

is composed of the Ti singlet ground state, the T\ ' triplet 
first excited states above the CF gap A = 0.7 meV and 
other excited states located at much larger energy above 

(2) 

~ 10 meV. The low-lying ri — T\ excitations should be 
responsible for the large 7 value. Furthermore, the study 
on Pr(Rui_a;Osa;)4Sbi2 series shows that a tiny variation 
of A leads to a large variation of 7 01 • 

In this Letter, we report on high resolution INS experi- 
ments on single crystals of PrOs4Sbi2, focusing attention 
on the low-lying singlet-triplet excitations. We found a 
clear softening of the excitation at the zone boundary 
Q = (1,0,0) and a peculiar Q-dependence of the excita- 
tion intensity. This result gives the direct evidence for 
the excitons derived from nonmagnetic quadrupolar in- 
teractions between 4f electrons. Furthermore, the tem- 
perature dependence of the excitations suggests that the 
new 'quadrupolar' excitons are coupled to the HF super- 
conductivity. 

Single crystals of PrOs4Sbi2 were grown by Sb-self- 
flux method. In order to gain counting statistics in the 
INS experiments, we used two single crystals. The to- 
tal mass is about 6.8 g. One of them is the same as 
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used for the previous INS experiments [10J . Their super- 
conducting transitions at T c ~ 1.8 K were confirmed by 
specific heat measurements performed on small pieces of 
each crystals. Neutron scattering experiments were car- 
ried out using the cold neutron three-axis spectrometers 
IN12 and IN14 located at the Institut Laue-Langevin, 
Grenoble. Data presented in this Letter were exclusively 
taken on IN14 and are consistent with IN12 prelimi- 
nary data. The incident beam was provided by a ver- 
tically focused pyrolytic graphite (PG) monochromator. 
A liquid-nitrogen cooled Be filter placed after the PG 
monochromator was employed to cut down higher-order 
contamination of neutrons. A horizontally focusing PG 
was used as analyzer. Energy scans at constant Q along 
the [1,0,0] direction were performed with the fixed final 
neutron energy 3.0 meV and with collimation open-60'- 
open-open. The energy resolution determined by the in- 
coherent signal is 0.10 meV full width at half maximum 
(FWHM). The assembly of two single crystals has a to- 
tal mosaicness of 1.3°, corresponding approximately to 
the mosaicness of each individual crystal. The assem- 
bly was mounted in a 3 He- 4 He dilution refrigerator and 
a standard 4 He cryostat with (hkO) as horizontal scatter- 
ing plane. The latter was used especially for measuring 
excitations in the vicinity of T c . A background was de- 
termined by using the neutron energy gain data and was 
subtracted. The intensity of neutron scattering along 
the [1,0,0] direction I bs(<2> is proportional to the dy- 
namic structure factor S(Q,u>) — S zz (Q 1 uj) + S yy (Q,uj) 
as l ohs (Q,uj) oc A{Q)f 2 (Q)S(Q,oj), where A(Q) is the 
absorption coefficient and f(Q) the magnetic form factor 
of Pr 3+ ion. A(Q) was determined by measuring the Q- 
dependence of the incoherent scattering intensity under 
the same geometrical condition for the scattering angle 
and the sample rotation angle as energy scans in Fig. ^ 
By using this experimentally determined absorption co- 
efficient, we got S(Q,u}). 

Figure n shows the dynamic structure factor S(Q,w) 
of PrOs4Sbi2 with different Q vectors along the [1, 0, 
0] direction at the lowest temperature 0.07 K and at 
3.9 K. The peaks around 0.4~0.8 meV correspond to the 

(2) 

Ti — CF transition. In both superconducting and 
normal states, the clear softening of the low-lying excita- 
tions at the zone boundary Q = (1, 0, 0) with minimum 
energies of 0.43 meV and 0.55 meV respectively has been 
observed. Note that the wave vector (1, 0, 0) is the same 
as that of the field-induced AFQ ordering, suggesting the 
importance of the quadrupolar degrees of freedom even 
in zero magnetic field, which is more directly evidenced 
by the (J-dependence of the excitation intensity as dis- 
cussed later. In the superconducting state, the spectra 
show sharp peaks with lincwidths of 0.2 meV, indicat- 
ing the localized character of 4f electrons. In the normal 
state, on the other hand, the peaks significantly broaden 
and their positions shift to higher energies. In Fig. 2] a 
faint peak seems to be around 1 meV. However, it has not 
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FIG. 1: (Color online) Energy spectra of excitations in 
PrOs4Sbi2 with different Q vectors along the [1, 0, 0] direc- 
tion at 0.07 K (open circle) and at 3.9 K (closed circle). Lines 
are the two Lorentzian fit. 



been observed in the IN14 experiment using the 4 He cryo- 
stat. Therefore, we consider that the slightly enhanced 
intensity at about 1 meV appears to be an artifact of the 
experimental set-up. The observed S(Q,uj) was fitted to 
two Lorentzian line shapes, corresponding to the tail of 
incoherent signal and the peak around 0.4^0.8 meV. For 
all the spectra, the fitting parameters of the incoherent 
tail are fixed as those determined at the lowest tempera- 
ture. 

Figure[21(a) is the dispersion curves obtained from the 
fitting in Fig.^ The energy of the excitations monoton- 
ically softens from the zone center to the zone boundary 
in both the normal and superconducting states, while the 
energy decreases uniformly for all Q vectors with decreas- 
ing temperature. This behavior seems to be accounted 
for as magnetic excitons, which have been extensively 
studied in rare earth metals and compounds for over 
three decades 0, but the interactions that modulate 
the CF levels to create so-called excitons are not necessar- 
ily magnetic exchange. Indeed, the agreement between 
the wave vector of the softening and the modulation vec- 
tor of the field-induced AFQ ordering |8| most naturally 
indicate that the interaction between quadrupolar mo- 
ments Qi remains the principal interaction. However, we 
must also take into account the interaction between mag- 
netic moments M.;, because the bulk magnetic suscepti- 
bility indicates finite antiferromagnetic correlations 0]. 
Therefore, the experimental determination of the nature 
of the main interaction is a fundamental issue to clarify 
the mechanism of the superconductivity in PrOs4Sb 12 . 

The monotonic decrease of the exciton energy from 
the zone center to the zone boundary suggests that the 
nearest- neighbor interaction is dominant in PrOs4Sbi2- 
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FIG. 2: (Color online) (a) Dispersion relations for the peaks 
around 0.4~0.8 meV and (b) Q-dependence of the integrated 
intensity along the [1, 0, 0] direction at 0.07 K (open circle) 
and 3.9 K (closed circle) obtained from S(Q,lu) in Fig. In 
Fig. |51(b), the intensity is normalized by that at Q = (2.2,0,0) 
for T — 0.07 K. Dashed lines are guides for the eyes 



We, therefore, assume the following effective Hamiltonian 



ofinteractions, H = D M £i,j M, • M 3 + Dq J2 i:j Qi ■ Qj 
where Dm and Dq denote the coupling constants 



|14j. wnere um 
of magnetic and quadrupolar interactions, respectively. 
The Fourier transforms of the coupling constants in bcc 
lattice are A(q) = 8A Hj =x , y , z cos(q J a/2) (i = M,Q). 
If magnetic and quadrupolar interactions coexist in the 
singlet-triplet system, the exciton dispersion relation 
within the random phase approximation is written as (lo| 



fajJr 



(A + 2|a| 2 J D M (q))(A + 2|/3| 2 J D Q (q)), (1) 



where a and (3 are the off-diagonal matrix elements of 

(2) 

magnetic and quadrupolar moments between Ti and I 1 ^ , 
respectively |l6(. When Dq = 0, Eq. becomes the 
usual form of magnetic excitons 13] . The observed dis- 
persion curves follow the antiferro-case of Eq. |Q, i-e. 
Dm,Dq > 0. As pointed out by Shiina et ai, the mag- 
netic and nonmagnetic interactions give the same energy 
dispersion but they give opposite Q-dependence of the 
intensity. Therefore, the intensity analysis is a definitive 
test that can distinguish the nature of the interaction. 
The theoretical integrated intensity I(q) = J* S(q, u>)dui 
is given by 0] 



I(q) oc 



/ A + 2|/?pj q (q) 
A + 2|a| 2 ^ M (q)' 



(2) 



When the antiferro-quadrupolar fluctuations give rise to 
the dispersion of excitations, Eq. predicts that I(q) at 
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FIG. 3: Energy spectra of excitations in PrOs4Sbi2 at 
Q — (1,0,0) for different temperatures. Lines are the two 
Lorentzian fit. 



the zone boundary, where the nonmagnetic quadrupolar 
fluctuations are stronger, is smaller than I(q) at the zone 
center. The opposite behavior comes from the fact that 
I(q) is proportional to the magnetic correlation function. 
This conclusion is universal and also does not depend on 
the details of system. 

In Fig. [5] (b), we plot the Q-dependence of the in- 
tegrated intensity determined by subtracting the inco- 
herent signal from S(Q,uj) in Fig. ^ an d summing the 
remaining counts. The observed Q-dependence of the 
intensity is quite different from that of the usual mag- 
netic excitons. The intensity at the zone boundary is 
smaller than that around the zone center in both normal 
and superconducting states. The experimental results 
agree with Eq. © for \f3\ 2 D Q > \a\ 2 D M . This result 
directly evidences that the observed excitonic behavior 
is derived mainly from nonmagnetic quadrupolar inter- 
actions. Within our knowledge, PrOs4Sbi2 is the first 
example of the excitons due to quadrupolar interactions. 

We also observed the strong temperature dependence 
of excitation spectra at Q — (1, 0, 0) as shown in Fig. 01 
By the same fitting procedure as in Fig. ^ we plot the 
temperature dependence of the energy and the linewidths 
(FWHM) of the peaks around 0.4-0.8 meV in Fig.H(a) 
and (b), respectively. Assuming that the interactions 
are only quadrupolar, the temperature dependence of 
the exciton energy is given by multiplying the term of 
Dq (q) of Eq. P by a temperature renormalization fac- 
tor zcxn ( -A/tflT ) ' m ^ e same wa y as usua l magnetic 
excitons [13. The theoretical curve for A = 0.7±0.1 meV 
and (3 2 Dq = 0.029 ±0.01 meV approximately reproduces 
the observed peak shift as shown by the dashed line in 
Fig. El (a). It is obscure whether the energy of the peaks 
decreases gradually or decreases with a distinct anomaly 
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FIG. 4: Temperature dependence of (a) the energy and (b) the 
linewidths of the peaks at Q — (1,0,0). Dashed line denotes 
the theoretical curve of the excitons (see text). 



at T c as decreasing temperature. This ambiguity may 
come from the fact that the superconducting energy gap 
is comparable with the CF gap A = 0.7 meV. On the 
other hand, the linewidths of the peaks exhibit a discon- 
tinuous decrease just below T c as shown in Fig. 0] (b), 
suggesting the close connection between the quadrupolar 
excitons and the superconductivity. This result agrees 
with the general observation that for rare earth ions dis- 
solved in superconducting materials, the CF linewidth 
reacts to the Cooper pairing 0, 0] . Exciton based su- 
perconductivity theory also predicts that the life time 
of excitons gets longer below T c [l^. The novelty of 
PrOs 4 Sb 12 is the exotic feedback between A and 7 due to 
the proximity of multipolar instability. The open under- 
lying questions are if the coherence among the excitons is 
associated with the entrance in a multipolar state of tiny 
order parameter and if there is any link of our neutron 
results with double structure of the specific heat at T c . 

The usual BCS superconductor PrRu4Sbi2 with nor- 
mal effective mass [12], |2Qj, |2l| also has the singlet-triplet 
CF level scheme but the CF gap between Fi and T\ 
5.26 meV, which has been determined by our recent INS 
experiments on a polycrystal sample 22] , is much larger 
than A of PrOs4Sbi2. This fact clearly indicates that 
the smallness of A, namely the existence of some de- 
grees of freedom in low energy region, is the necessary 
condition for the mass enhancement. This is consistent 
with a mechanism of the mass enhancement through the 
c — / magnetic exchange in the singlet ground state sys- 
tem [2^. In addition, when 4f electrons have magnetic 
degrees of freedom in such low energy region, the local 
magnetic fluctuations may be expected at low tempera- 
tures. In such a case, the resulting magnetic quasi-elastic 
scattering is usually detected in INS experiments as in 
Ce-based HF compounds. In PrOs4Sbi2, however, quasi- 
elastic scattering has not been detected at least above 
0.2 meV as shown in Figs.n an d[3] This result suggests 
other nonmagnetic mechanism rather than that due to 



magnetic fluctuations. On the other hand, the significant 
broadening of spectral shape in the normal state clearly 
suggests the importance of the c — / interaction. Com- 
bining these facts, the mass enhancement as well as the 
superconductivity in PrOs4Sbi2 may originate from the 
strong c — / interaction through the quadrupolar chan- 
nels. 

In summary, new 'quadrupolar' excitons, which are CF 
excitations modulated by quadrupolar interactions, are 
realized in PrOs4Sbi2- A clear softening of this excita- 
tion occur at Q = (1,0,0). On cooling, all the exciton 
branch softens. Moreover the tendency of the narrowing 
of linewidths below T c suggests that the quadrupolar ex- 
citons are related to the HF superconductivity and vice 
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